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The potential for a precise measurement of the Unitarity Triangle angle 7 in future experiments 
from the decay B° DK*'^ is well-known. It has recently been suggested that the sensitivity can 
be significantly enhanced by analysing the -B" — > DK'^-k~ Dalitz plot to extract amplitudes relative 
to those of the flavour-specific decay K^ . An extension to this method which includes 

the case where the neutral D meson is reconstructed in suppressed final states is presented. The 
sensitivity to 7 is estimated using this method and compared to that obtained using the B° DK*° 
decay alone. Experimental effects, such as background contamination, are also considered. This 
approach appears to be a highly attractive addition to the family of methods that can be used to 
determine 7. 
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I. INTRODUCTION 

The quark flavour sector of the Standard Model of 
particle physics, described by the Cabibbo-Kobayashi- 
Maskawa (CKM) quark mixing matrix [TJ gives a 
successful description of all current experimental mea- 
surements of quark flavour-changing interactions. It 
also provides an excellent laboratory to search for ef- 
fects of physics beyond the Standard Model (see, for 
example, Refs. [S] HI |5]). A critical element of this 
programme is the precise measurement of the an- 
gle 7 = arg i-VudV*^/VcdV*^) of the Unitarity Triangle 
formed from elements of the CKM matrix. 

A method to measure 7 with negligible theoretical un- 
certainty was proposed by Gronau, London and Wyler 
(GLW) jSHZl. The original method uses B DK decays, 
with the neutral D meson reconstructed in CP eigen- 
states. The method can be extended to use D meson 
decays to any final state that is accessible to both 
and D^, in particular doubly-Cabibbo-suppressed decays 
such as if+TT" lUEl have been noted to provide enhanced 
sensitivity to CP- violation effects. 

The use of neutral B decays is particularly interesting 
since the two contributing amplitudes are more similar in 
magnitude, so that direct CP- violation effects may be en- 
hanced relative to those in charged B decays. The decay 
P" DK*'^ is especially advantageous since the charge 
of the kaon in the K*'^ K^Tr~' decay unambiguously 
tags the flavour of the decaying B meson, obviating the 
need for time-dependent analysis |10j . This appears to 
be one of the most promising channels for LHCb to make 
a precise measurement of 7 [TTJ [T5J [T3] . 

The approach that has mainly been considered until 
recently is a quasi-two-body analysis of P° DK*^ . In 
this analysis, the contributions from other resonances in 
the P° DK^TT^ Dalitz plot that interfere with the 
K*''^ within the selected mass window are handled by the 
introduction of an additional hadronic parameter |14j . 



This parameter, normally denoted by k, takes values be- 
tween and 1 where implies that all sensitivity to 7 
is lost, and the limit of 1 is reached in the case that 
no amplitudes other than DK*'^ contribute. Estimates 
suggest that 0.9 < k < 1.0 for a K*'^ mass window of 
±50 MeV dS]. 



Recently it has been noted that the natural width 
of the K* meson can be used to enhance the sensitiv- 
ity to the CP- violating phase 7 through analysis of the 
B^ DK^TT^ Dalitz plots [T5J [T7]. By comparison of 
the Dalitz-plot distributions of events in the cases where 
the neutral D meson is reconstructed in flavour-specific 
and CP-eigenstate modes, the complex amplitudes of 
the DK*^ decays can each be determined relative to the 
flavour-specific amplitude. This allows for a di- 

rect extraction of 7 from the difference in amplitudes, 
rather than from the rates. 



In this paper we extend the method proposed in 
Ref. [17] to include also the case where the neutral D 
meson is reconstructed in suppressed final states. This 
allows us to make a direct comparison of the sensitivity 
to 7 between the quasi-two-body analysis and the Dalitz- 
plot analysis. We also study possible systematic effects 
that may limit the sensitivity of the analysis, including 
uncertainties on the correct composition of the Dalitz- 
plot model and a brief discussion of experimental effects. 



The remainder of the paper is organised as follows: in 
Section [ll] we give an overview of the method; in Sec- 
tion |III| we describe the Dalitz-plot model used for our 
study; in Sections |IV| and |V] we present the results we 
obtain in the quasi-two-body analysis and in the Dalitz- 



plot analysis, respectively; in Section VI we discuss how 
experimental effects can be handled, before summarising 
in Section rVlTl 
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II. METHOD 

In this Section we describe the various methods that 
can be used to extract 7 from B —> DK'^*\ We first 
review the quasi-two-body approach, and then recap the 
recently proposed Dahtz-plot technique [17,1 before de- 
scribing the extension to the method used in our study. 
All methods exploit the interference between the two 
tree-level amplitudes shown in Fig. [l] Conventionally, 
the ratio of magnitudes of these two amplitudes is re- 
ferred to as tb, while their strong phase difference is 
labelled 5b- 



w 




B" 



W 




FIG. 1: Feynman diagrams for DK*° , via (left) a 

h — > cus transition and (right) a 6 — > ucs transition. 

The GLW method [H of extracting 7 uses the fol- 
lowing rates and asymmetries in i? ^ DK^*^ decays: 



R+ = 



A+ = 



T(B ^ Di^^Ki*)) + T{B ^ Df,,K(*)) 
1 + r| ± 2 cos((5b) cos(7) , 
T(B D±K^*^) - T{B D±K^*^) 
T{B D±K(*)) + T{B D±K(*)) 
±2rB sm{SB) sin(7) 



(1) 



(2) 



Here, B (B) is used to refer to either B or or 
B'^), while D± refers to a neutral D meson reconstructed 
in a CF-even {eg. K+R-) or CP-odd (eg. if^TrO) final 
state, and Z?fav refers to a neutral D meson reconstructed 
in a favoured, quasi-flavour-specific {eg. D'^ —f K^tt^) 
final state. Note that experimentally it is convenient 
to measure R± normalised to an equivalent double ra- 
tio from B Dtt or B ^ Dp decays. 

Since R+A^ + R_A_ ~ 0, the above four observables 
give three independent constraints on the three parame- 
ters 7, rB and 5b. This is sufficient to solve the system up 
to an eightfold ambiguity. However, when measurements 
are performed in a hadronic environment as at LHCb, 
the reconstruction of the CP-odd final states becomes a 
significant experimental problem. Therefore, additional 
observables are required. 

A solution is to include doubly-Cabibbo-suppressed D 
meson decays, as first suggested by Atwood, Dunietz and 
Soni (ADS) |Hlig. Due to the fact that the D meson de- 
cays to K^TT^ are not truly flavour-specific, but include 
a suppressed contribution which is given by tdR^^" rel- 



ative to the favoured decay amplitude,^ enhanced CP- 
violation effects can occur. The rates and asymmetries 
of the B DK^*") decays to the suppressed final states 
are then given by 



Pads — 



r(P DsupK(*^) + r{B ^ D,^pK(*)) 

T{B ^ Aavif + r(P ^ Aavif 



= r% + rl, + 2rBrD cos{5b - 5d) cos{j) , (3) 



A 



ADS — 



r(p 

2 rB ro s\n{5B 



supK^*^) - T{B 
supK(*'>) + T{B 



D. 



5d) sin(7) 



Pads 



(4) 



Since the hadronic parameters of the D decay (r^j and 
5d) can be determined independently [TSl[Tn], these mea- 
surements provide two additional linearly independent 
constraints that can be used in combination with the 
GLW observables to obtain bounds on the three un- 
known parameters 7, r_B and 5b. With the decay modes 
D K^TT^, the ADS observables are well suited to re- 
construction in a hadronic environment. Consequently, 
one of the most promising strategies for the tree-level 
determination of 7 at LHCb is that from the combina- 
tion of measurements of P+, A+, Pads and Aads in 
charged B DK decays or neutral B DK* de- 
cays nap. 

The finite width of the K*'^{892) resonance leads to 
additional complications in the analysis of P ^ DK* de- 
cays, since other contributions to the B DKtt Dalitz 
plot can affect the population within the K* mass win- 
dow. This can be handled by making the following sub- 
stitutions [Ml: 



rB — > rs 



/dp I^mI dx 



\/dp 



\Aydx 



/pp|A„||A,|e'^df 
/dp \Au\^ dxj^^ \Ac\^ dx 



(5) 
(6) 



where Ac and A^ are respectively the amplitudes car- 
rying the phase of the b cus (Fig. [l|left)) and of the 
h ucs (Fig. [ijright)) transitions, and 6 is the strong 
phase difference between them, all as functions of the 
Dalitz-plot position x. The integrals are over the region 
of the Dalitz plot that is defined as the K* mass win- 
dow. In the limit that DK* is the only contribution in 
this window, rs — > rB, 5s — > 5b and n — > 1. 

With this treatment, the two hadronic parameters as- 
sociated with the DK* decay (r^ and 5b) are replaced 
with two effective parameters (rg and 5s) and a new un- 
known (k) is introduced. Since the combination of (CP- 
even) GLW and ADS observables provides four linearly 



^ Note that the sign convention for 60 used in this paper is op- 
posite to that used in most of the literature on 7 measurements. 
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independent measurements, it is possible to determine k, 
from the data together with 7, and 5b- Alternatively, 
external input, either theoretical or experimental, could 
be used to constrain n [15] , 

We refer to the extraction of 7 using the approach out- 
lined above as the quasi- two-body analysis. The addition 
of ADS observables helps to resolve two of the ambigui- 
ties of the GLW approach [21^; however, the effectiveness 
of this depends on the values of ru, 5d and Jb, as well as 
the statistical sensitivity. Moreover, the sensitivity to 7 
depends on the values of the unknown hadronic parame- 
ters, particularly 5b- 

The recently proposed B DKtt Dalitz-plot analy- 
sis [17] exploits the presence of the B contri- 
bution that serves as a reference amplitude, since the 
flavour of the neutral D meson produced in — > Dtt^ 
is tagged by the charge of the accompanying pion. Con- 
sidering flavour-specific D mesons, we can define the 
B DK* amplitude relative to this reference, as il- 
lustrated in Fig.|2];ieft), 



^(go ^ d0k*°) 

A{B" d;-k+) 



ge 



(7) 



Note that in Eq. [7] and throughout the discussion 
below we neglect factors of A{K*'^ K^Tr~) and 
A{D2^ D^TT^) that formally should appear in the nu- 
merator and denominator respectively, since they even- 
tually cancel in the observables of interest. 

Considering now CP-even D mesons, using the con- 
D° ±W), we find (see Fig. [^Fright)) 



vention D+ = -h= 

v2 



^/2A{B^ ^ D+K*°) 
V2A{B^ 



d;-k+) 



ge 



1 + tbc 



i(i5B+7) 



) , (8) 



where £'2+ denotes that the neutral D meson produced 
in the decay of the I?2~ is reconstructed in a CP-even 
eigenstate. Thus, we find [T7] 



{V2A{D+K*^))/{V2A{Dt 



K+)) 



(A{DOK*^))/{A{D*-K+)) 
^A{D+K*^) 



A{D°K*^) 



(9) 



where the variables (a;+,y+) are the same as those used 
in the analysis of P — > DK with D — 
cays 1221 [23]. Constraints on + iy^ 



tbc 



is. 



de- 

i-7) 



are likewise obtained from equivalent expressions for the 
charge-conjugate P° decays. The extraction of 7 from 
this Dalitz-plot analysis with only a single unresolvable 
ambiguity (7 — > 7 -I- tt, i5b — > 5b + t^) is possible using 
only CP-even D decays. Consequently, we restrict our 
discussion to CP-even D decays, since those are experi- 
mentally accessible in a hadronic environment; however, 
we note that CP-odd decays give similar expressions, but 
with the right-hand side of the last two relations of Eq. [9] 
multiplied by a minus sign. 



The discussion above, and in Ref. jT7|, considers that 
the D mesons used for normalisation are reconstructed 
in flavour-specific final states. However, as already men- 
tioned, the favoured decay D'^ —> K~Tr~^ is only approx- 
imately flavour-specific. Furthermore, since the doubly- 
Cabibbo-suppressed decays are used to great benefit in 
the quasi- two-body analysis, it is reasonable to ask if 
they can also be included in the Dalitz-plot analysis. We 
therefore extend the method to include the effects of the 
suppressed I3-decay amplitudes. We find 



A(P° ^ AavA-*P) 

A{B° ^ Aupi^*") 



= ge 



yl(P" 



(10) 
(11) 



while the expression for the amplitude ratio in the Dalitz 
plot with the CP-even D meson is unchanged from Eq. |8] 
These amplitudes are illustrated in Fig. [3] As before, 
the equivalent expressions for the charge-conjugate B'^ 
decays are obtained with the substitution 7 — > —7. 

As one would expect, the amplitudes obtained from 
the Dfav-fsT+TT^ Dalitz plot (Fig. 3]^left)) are hardly dis- 
tinguishable from those for the idealised flavour-specific 
D decays. As discussed in Ref. [T^, if the suppressed 
amplitudes are neglected, it will lead to only a small bias 
in the extraction of 7. 



The amplitudes obtained from the -DsupA 



Dalitz 



plot (Fig. [3]^ right)) are markedly different from those of 
the other Dalitz plots, and this new information can in 
principle be used to constrain 7. However, in this case 
the D2K contribution is no longer suitable as a reference 
amplitude, due to its small size, and the measurement 
of the relative phase between D2K and DK* amplitudes 
would be expected to have a large uncertainty. Nonethe- 
less, it should be possible to obtain information about 
the relative magnitude of these amplitudes, which will 
provide sensitivity to 7. Note that, in contrast to the 
quasi- two-body analysis, it is not necessary to use exter- 
nal constraints on ru and 5d in the Dalitz-plot analysis. 



III. 



DKtt DALITZ PLOT MODEL 



We construct a model for B DKtt Dalitz plot dis- 
tributions using the isobar formalism, in which the total 
amplitude is written as the coherent sum of contributions 
from resonant and nonresonant terms: 

M{x) = anre*'^"-f (12) 

-^Rb{x)FR^did2ix)BWr{x)Sr{x)- 



In Eq. 12 x = (™/f^,w|,^) represents the position in 
the Dalitz plot, ae'^'^ describes the complex amplitude for 
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FIG. 2: Argand diagrams illustrating the measurements of relative amplitudes and phases from analysis of the Dalitz plots of 
(left) D°K^Tv~ and (right) D+K^tv~ . In these illustrative examples the following values are used: g = 1.5, A — 20°, 7 = 75°, 
5b = 45° and rs = 0.4. 



ImA 




IrriA 



-1 



e-'^"A(B''^D3„„K°) 




■-\e-''"A(B<'^D;-3,pn 



FIG. 3: Argand diagrams illustrating the measurements of relative amplitudes and phases from analysis of the Dalitz plots of 
(left) Dia.vK^ -K^ and (right) Dsupif^Tr". Note that in the latter the amplitudes are rotated by —5d to maintain the convention 
of having the DJ~A'"'" amplitude on the real axis. In these illustrative examples the following values are used: g — 1.5, A = 20°, 
7 = 75°, 5b = 45°, TB = 0.4, 5d = -158° and m = 0.06. 



each component, the F terms denote vertex form factors, 
BW the resonance propagator and S the Lorentz invari- 
ant spin factor. We use Blatt-Weisskopf barrier form 
factors [24, and use relativistic Breit-Wigner hneshapes 
to describe the propagators. We use the Zemach formal- 
ism [25l [26] for the spin factors. We assume that the 
nonresonant contribution is constant across the phase 
space. This is a sufficient approximation for the study 
at hand, even though a more complicated description is 
hkely to be necessary to fit real data. All amplitudes are 
evaluated using the qft++ package [27]. 

We develop a model for the DK^tt^ Dalitz plot 



distribution based on the following results from Ref. [55] : 

BiB° ^ DK+TT-) = (88± 15±9) X 10-^ 
B{B" ^ DK*°{892)[K+TT-]) = (38 ± 6 ± 4) x 10"^ 
B{B° ^ D;-{2460)[DTr-]K+) = 

(18.3 ±4.0 ±3.1) X 10-^ 
B{B^ ^ DK+TT")^, = (26±8±4) X lO"*'. 

The results for the resonant contributions were extracted 
using events in the regions \Mk-k — Mk-'« \ < 150 MeV and 
l^fcir — -^^D'- I < 75 MeV, respectively. The strengths 
of the if* (892) and 152" (2460) ampUtudes in our model 
were set by requiring that the fit fractions of these reso- 
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TABLE I: Parameters of resonances used in the model |30| . 



Resonance 




Mass (MeV) Width (MeV) 


A'* (892) 


1" 


896 


51 


K5{U30) 


0+ 


1412 


294 


7^2* (1430) 


2+ 


1432 


109 


74:*(1680) 


1" 


1717 


322 


L>S(2400) 


0+ 


2403 


283 


D5(2460) 


2+ 


2459 


25 


L>:i(2700) 


1" 


2690 


110 



nances in these regions match the pubhshed results [25] . 

Since we expect the Dahtz plot to contain contribu- 
tions from other resonances not considered in Ref. (25] . 
we add additional contributions based on results from 
an analysis of the B° D^K^t:^ Dalitz plot [29], tak- 
ing isospin and colour-suppression factors into account 
as appropriate. Some additional scaling of these reso- 
nance terms was performed to provide a better match to 
the results published in Ref. [28]. These experimental re- 
sults provide information about the likely magnitude of 
contributions from K* and D* resonances to the Dalitz 
plot. 

Additionally, contributions from £'*-type resonances 
can in principle contribute. Their effect could be sig- 
nificant since they are mediated by 6 — > m transitions 
which provide the sensitivity to 7. The Z?*2(2573) and 
£'*]^(2700) states are known to decay to DK, and the 
latter has been observed in B decays [3T]. These are not 
included in our nominal model, but we consider their po- 
tential effect among the model variations discussed below 
in Section [V] 

Table [l] summarises the parameters of the resonances 
used in our analysis. Table [H] gives the fit fractions and 
relative phases of the various intermediate states that 
contribute to our B DKtt Dalitz plot model. The pa- 
rameters given are appropriate for the case that the neu- 
tral D meson decays to a favoured, quasi-flavour-specific 
final state (namely, DK^ti^ , D — > K^tt' and 

charge conjugate). An example Dalitz-plot distribution 
generated from this model is shown in Fig. [4] 

We derive the Dalitz-plot distributions for other D- 
decay final states using the equations presented in Sec- 
tion |ll] and initially taking 7 = 60°. We consider several 
different possible values of tb and 5b for the DK*^{^%2) 
amplitudes, while for the less significant contributions 
from other DK* channels we simply set rs — 0.4 and 
5b = 0°. We note that our nominal value of = 
0.4 is at the upper end of the experimentally allowed 
range for this parameter (e.g. the UTfit Collaboration 
gives rB{DK*"{892)) e [0.081,0.397] at 95% confidence 
level [32 ), but it is convenient to use this value to allow 
comparison with previous studies. 

We use our Dalitz-plot model to generate ensembles 
of event samples corresponding to the different D meson 
final states. The numbers of events are based on ex- 
pectations for one year of nominal luminosity at LHCb 
(2 fb~^) [121, assuming that the efficiency does not vary 



^ 25 




2 4 6 8 10 12 

m|,(GeV^) 

FIG. 4: High-statistics B Di^vKn Dahtz plot distribu- 
tion generated with our nominal model. Structures due to 
DK*°{892), DKfiUSO) and D^' {2-i60)K+ , as well as the 
nonresonant contribution, are clearly apparent. 

TABLE II: Summary of the B Dt^^Kix Dalitz plot model. 



Intermediate state Fit fraction Phase (°) 



D7^*"(892) 


0.46 







0.0001 


284 




0.12 


221 


7? A'*" (1680) 


0.02 


128 


D;-(2460)A'+ 


0.34 


325 




0.02 


267 


nonresonant DKn 


0.06 


140 



across the Dalitz plot. The exact numbers in each sample 
vary as functions of the input values of the parameters 
(particularly 5b and 7), but typically are around 7300 
for Aavi^TT, 700 for D^upKir and 600 for D+Kn (for B° 
and B° decays combined). 

IV. RESULTS WITH THE QUASI- TWO-BODY 
APPROACH 

We first study the sensitivity to 7 in the quasi-two- 
body approach. We take the six Dalitz-plot distributions 
(DfavivTTT, DsupKTT and D+Ktt all for both W and B°) 
generated as described in the previous section, and apply 
selection requirements to select the DK* dominated re- 
gion. We then perform minimisation to fit the yields 
in each sample to determine five parameters: j,rs,5s, k 
and an overall normalisation. The parameters td and 5d 
are fixed to their measured values (0.0616 and —158°, re- 
spectively). For each pseudo-experiment we perform ten 
fits with initial values of the parameters randomised in 
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the ranges 7 e [0, 2tt), tb € [0, 1), 6b G [0, 27r), k G [0, 1), 
and we take the results of the fit with the smallest ■ 

The selection requirements for the K* are an interest- 
ing sub-topic worthy of some discussion. The width of 
the Kit invariant-mass window around the nominal K* 
mass affects the size of the event samples. Increasing the 
width yields greater statistics so that one would naively 
expect a reduction in the statistical uncertainty on 7. 
Unfortunately, increasing the width of the Ktt invariant- 
mass window also leads to an increase in the dilution 
from the non-DK* component of the amplitude, i.e. a 
decrease in k, which tends to decrease the sensitivity to 
7. Hence, it is important to optimise the K* selection 
requirements. 

One can increase k by introducing a £'2(2460) 
veto (e.g. rejecting events that satisfy 

\Md-k - Md*\ < 75 MeY). This veto has a rather 
minimal impact on the statistics. In our nom- 
inal model with Sb — 180°, choosing to use 
- Mk'{892)\ < 150 MeV yields k = 0.93. 
Applying the £'2(2460) veto above increases this to 
K — 0.97. Decreasing the Ktt invariant-mass window to 
\M. 



Ktt 



if* (892) I 



< 50 MeV yields k — 0.99; however, 
with the limited statistics in our pseudo-experiments, 
this decrease in dilution is counter-balanced by the loss 
of statistics. Experiments with higher statistics may 
benefit by using a tighter Ktt invariant-mass window. In 
the results below, we proceed using the -D2(2460) veto 
discussed above and requiring the Kn invariant mass to 
be within 150 MeV of the nominal K* mass. 

Results from the pseudo-experiments for various dif- 
ferent input values of Sb are shown in Fig.|5] It is appar- 
ent that there is a strong correlation between the fitted 
values of 7 and ^5, and moreover that there are ambi- 
guities in the solution that are particularly pronounced 
for values of Sb near 90°. The locations of each of the 
ambiguous solutions for the case where ru — are also 
shown in Fig. |5] The ambiguities in the solutions found 
by fitting our data are close to these values. Since r^i is 
small (our data were generated with ru = 0.0616), sta- 
tistical fluctuations can lead to the best value being 
obtained near one of the ambiguous solutions. The in- 
ability to resolve such ambiguities is a limitation of the 
quasi- two-body method. 

As a consequence of these ambiguities, the distribu- 
tions of the fitted values of 7 are not Gaussian. There- 
fore, in Tab. |III| we report both the mean and rms of 
each distribution as well as the corresponding values ob- 
tained by fitting the data in the regions near the gener- 
ated values of 7 to Gaussian lineshapes, i.e. ignoring the 
ambiguities in the solution. For 6b = 90°, it is not possi- 
ble to separate the correct solutions from the ambiguous 
solution, so we have not performed a Gaussian fit to the 
distributions for this sample. 

We have also considered an alternative approach to the 
fits, in which the value of k is fixed. Although n can only 
be measured from analysis of B ^ DKtt Dalitz plots (so 
that if this measurement can be performed, the Dalitz 



plot analysis will most likely be feasible), it is conceiv- 
able that K may be determined from theory, or from a 
different experiment. We therefore repeat the fits with 
K constrained to the values obtained directly from our 
model. This approach leads to an improvement in the 
resolution on 7 near some solutions {e.g., near the true 
solution for Sb — 0°), but does not remove the ambi- 
guities (see Fig. |6]). The results obtained by fitting the 
regions near the generated values of 7 to Gaussian line- 
shapes can be found in Tab. |III[ This provides a compar- 
ison with Ref. [12], where the same approach was used 
to estimate the sensitivity to 7 using the quasi-two-body 
approach on one nominal year's data from LHCb. De- 
spite several important differences between our toy study 
and the detailed sensitivity study of Ref. [12 , the results 



are in very good agreement (see Tab. Ill I. This provides 
confidence in the absolute scale of our sensitivity esti- 
mates. 



TABLE III: Results for 7 obtained from the quasi-two-body 
analysis, with different input values oi Sb and a comparison 
with the results of Ref. [12]. These results are obtained with 
7 = 60° and rs = 0.4. See text for details. 







n free 


7(° 


K locked 


Ref. [H] 






Distribution 


Gaussian fit 


Gaussian fit 






mean rms 




a 




O" 


a 


Sb = 


= 0° 


55.5 14.6 


57.7 


7.3 


60.1 


5.9 


6.2 


Sb = 


45° 


60.2 19.6 


57.4 


10.6 


57.8 


10.3 


10.8 


Sb = 


90° 


88.3 17.9 










12.7 


Sb = 


135° 


60.4 18.0 


56.6 


9.4 


56.3 


9.7 


9.5 


Sb = 


180° 


55.1 19.2 


56.6 


7.4 


59.7 


5.5 


5.2 



The results obtained by varying in our model are 
given in Tab. jlVj As is well known, the sensitivity to 7 
reduces as becomes smaller. It is also clear that the 
effect of ambiguities becomes more significant, since the 
discrepancy between the rms of the whole distributions 
and the width of the Gaussian peak near the correct so- 
lution increases. 



TABLE IV: Results for 7 obtained from the quasi-two-body 
analysis, with different input values of rs. These results are 
obtained with 7 = 60°, Sb = 0° and k free. 







in 






Distribution 


Gaussian fit 






mean rms 


/I a 


TB 


= 0.4 


55.5 14.6 


57.7 7.3 


TB 


= 0.3 


57.0 21.1 


56.8 9.7 


rB 


= 0.2 


58.9 29.3 


56.9 11.8 



V. RESULTS WITH THE AMPLITUDE 
ANALYSIS 

To study the sensitivity to 7 using the Dalitz-plot anal- 
ysis, we perform a simultaneous maximum likelihood fit 
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FIG. 5: Distributions of fitted values of 7 in the quasi-two-body approach, with different input values oi 5b- These results are 
obtained with 7 = 60° and tb = 0.4. The location of each of the generated solutions for k= \, i.e. 5s ~ 5b, is given by a red 
solid cross. The locations of each of tlie ambiguous solutions obtained for k = 1 and ro ~ are given by the intersections of 
the blue dashed Hues. 
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FIG. 6: Distributions of fitted values of 7 in the quasi-two-body approach with k free (left) and fixed (right). These results 
are obtained with 7 = 60°, 5b = 90° and rB = 0.4. 



to all six Dalitz-plot distributions. We initially fix rjj and 
Sd to their measured values, leaving 24 free parameters 
(7, four parameters - A, g, rs ^ Sb - for each ii'*-type 
resonance and nonresonant contribution, and two ~ mag- 
nitude and phase - for each Z?*-type resonance, with one 
phase fixed). As before, for each pseudo-experiment we 
perform multiple fits with randomised initial parameter 
values, and we take the results of the fit with the smallest 



negative log-likelihood. 

In Section |IV| we found that the quasi-two-body ap- 
proach was, in many cases, unable to resolve ambiguities 
in the solution. This was especially true for values of 
5b near 90° . Figure [7] shows the distributions of results 
from the pseudo-experiments generated with 5b — 90° 
for the amplitude and quasi-two-body approaches. The 
additional information utilised in the amplitude analysis 
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FIG. 7: Distributions of fitted values of 7 in the amplitude (left) and quasi-two-body (right) approaches. These results are 
obtained with 7 — 60°, Sb = 90° and rs = 0.4. 



is sufficient to remove the ambiguities in the solution. 

The distributions of the fit results for the parameters 
rB{DK*{892)), Sb{DK*{892)) and 7 obtained from the 
pseudo-experiments generated with Sb = 0° are shown in 
Fig. [8] We fit the distributions with Gaussian lineshapes, 
since all are consistent with this shape (this is true for all 
Sb values), and report the means and widths in Tab. [V] 
The resolution on 7 varies between 4.2° — 5.8°, depending 
on the value of , in our nominal model. Thus, unlike for 



the quasi-two-body approach (Tab. Ill), the sensitivity to 
7 is not strongly dependent on the value of Sb ■ 



TABLE V: Results for r b{DK* (892)), Sb{DK* (892)) and 7 
from the Dalitz-plot analysis, with different input values of 
Sb- These results are obtained with 7 = 60° and rs = 0.4. 







rB 
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/i a 


Sb = 


= 0° 


0.40 


0.01 


0.8 4.9 


59.8 4.4 


5b = 


45° 


0.40 
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46.2 6.0 


59.2 5.5 


Sb = 


90° 


0.40 


0.01 


90.2 6.4 


59.8 5.7 


Sb = 


135° 


0.40 


0.01 


134.0 6.3 


59.3 5.8 


Sb = 


180° 


0.40 


0.01 


179.8 4.2 


59.7 4.2 



The physical values of and 7 may be different from 
those used in our model. To study how this would af- 
fect our results, we vary the values of and 7 in our 
model and fit the data using the same procedure outlined 
above. The results obtained considering 0.2 < < 0.4 
are given in Tab.|VI| As expected, the sensitivity to 7 de- 
creases as tb decreases; however, even for tb = 0.2 the 
distribution of fit results is still approximately Gaussian, 
and the resolution on 7 is still 7°. The results obtained 
considering 45° < 7 < 75° are given in Tab. | VII] We con- 
clude that the sensitivity to 7 does not depend strongly 
on the value of 7 itself. 

There are a number of parameters in our model which 
are not well constrained by data; however, most of these 
are not expected to impact significantly the sensitiv- 
ity to 7. For example, the strong phase between the 
D;-(2460)K+ and Z)X*0(892) amplitudes, which we de- 
note by A [D if*" (892)], could be very different from the 



TABLE VI: Results for tb^DK* {892)), 5s (DA'* (892)) and 
7 from the Dalitz-plot analysis, with different input values of 
Tb- These results are obtained with 7 = 60° and Sb =0°. 
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TABLE VII: Results for TBiDK* {892)), (5b (DA"* (892)) and 
7 from the Dalitz-plot analysis, with different input values of 
7. These results are obtained with tb ~ 0.4 and Sb ~ 0° . 
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0.40 0.01 


0.4 4.5 


74.8 4.0 



value used in our model. To determine what effect the 
value of this parameter has on the extracted value of 7, 
we vary A[DK*'^{892)] in our model over the physically 
allowed range, [0, 2tt), and refit the data. The fit results, 
given in Tab. |VIII| show that the impact on the sensitiv- 
ity to 7 is minimal. 



TABLE VIII: Results for rs(DA'*(892)), 5b{DK*{892)) and 
7 from the Dalitz-plot analysis, with different input values of 
A[£)A'*°(892)]. These results are obtained with rs = 0.4, 
Sb = 0° and 7 = 60°. 
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0.40 0.01 


0.5 4.9 


59.7 4.1 


145° 


0.40 0.01 


0.4 4.9 


59.7 4.3 


235° 


0.40 0.01 


0.0 5.0 


59.5 4.0 


325° 


0.40 0.01 


0.8 4.9 


59.8 4.4 



None of 



for the 



(1430), 



the parameters 
DX*°(1430) or DK*^{\&8IS) amplitudes is well con- 
strained by data. To see what effect the presence of CP 
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FIG. 8: Distributions of fitted values of rs (left), 5b (middle), and 7 (right) in the amplitude analysis. These results are 
obtained with 7 — 60°, Sb = 0° and rs ~ 0.4. 



violation in these amplitudes has on the sensitivity to 7, 
we set rB{DK^°{U30)) = 0, rB{DKf{U30)) = and 
rB{DK*° {1680)) = in our model and refit the data. 
The results obtained for this model variation are given 



in Tab. IX There is no appreciable bias on the extracted 
value of 7, and the resolution on 7 decreases by only 
1.2°. This is not surprising given the relatively small 
contributions to the total B —> DKtt amplitude from 
DK^°{U30), DK^°{U30) and DK*°{1680) (see Tab.|ll|. 
If, however, in real data these or any other DK*-type am- 
plitudes do contribute strongly to i? — > DKtt (and have 
sufficiently large CP violation), then the sensitivity to 
7 could be better than that obtained from our nominal 
model. To further test how mismodelling of the Dalitz 
plot may affect the extracted value of 7, we refit the 
data from our nominal model using only the DK* (892), 
D^(2460)iir and nonresonant amplitudes (fixing all oth- 
ers to zero). The results are shown in Tab. |X] Again, 
there is little impact on the extracted values, indicating 
that this analysis may suffer much less model-related un- 
certainties than some other methods to extract 7 |22l ES] . 



TABLE IX: Results for rfl(D_ft:*(892)), (5s (DA" (892)) and 
7 from the Dalitz-plot analysis, with r_g(-DA"o''(1430)) = 0, 
rB{DKf(U30)) = and r b{DK*" (1680)) = 0. These re- 
sults are obtained with rs — 0.4, Sb ~ 0° and 7 = 60°. 
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7(°) 
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H a 


0.40 0.01 
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5.2 


59.0 5.6 



TABLE X: Results for rs(DJ<'*(892)), (5b(DA'*(892)) and 
7 from the Dalitz-plot analysis, where the fit model con- 
tains only the Di('*(892), _DJ(2460)A' and nonresonant am- 
plitudes. These results are obtained with tb ~ 0.4, 5b = 0° 
and 7 = 60°. 
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0.40 0.01 
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We have also tested how the addition of /?s-type res- 
onances may affect the analysis. We add to the model a 
contribution from Z?*j^(2700), the magnitude of which is 
based on the assumption 



B{B 



■ Dsjtt) 
DsjD) 



B{B 
B{B 



■ Di*^n) 
Di*^D) 



(13) 



The right-hand side of this expression is found to be 
(2.4±0.6) X 10-3 for and (3.3±0.9) x IQ-^ for D* 
Using results from Ref . [31] , we estimate 

BiB -> i:>*i(2700)7r) x B{D;^{2700) ^ DK) w (14) 
2.7 X 10^3 X 11.3 X 10""* w 3 X 10"*^ . 

We add a contribution from this amplitude, with an 
arbitrary phase, to the model described in Tab. |TT] The 
results of fits to pseudo-experiments generated with this 
modified model are summarised in Tab. IXII There is lit- 
tle or no improvement in the sensitivity to 7 compared to 
the nominal model. This can be understood since the cor- 
ner of the Dalitz plot where the D*i(2700) and K*{892) 
would interfere is excluded by kinematic constraints (see 
Fig. [4]). Following this reasoning, larger effects from D^j- 
type resonances would be possible if there were significant 
contributions to the Dalitz plot from heavier K* and/or 
heavier D^j resonances. 



TABLE XI: Results for rs(L»iC* (892)), 5b{DK*{892)) and 7 
obtained from the Dalitz-plot analysis, where the Z)*i(2700) 
has been added to our nominal model. These results are for 



rs = 0.4, 5b = 0° and 7 = 


60° 
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7(°) 
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0.40 0.01 


0.8 


5.0 


59.9 4.3 



In our nominal fit, we fix the values of and ru. How- 
ever, the value of is currently not precisely measured: 
5^, = (-158±11)° [18lll9l|33]. To study what effect this 
parameter can have on our results, we rerun the fits de- 
scribed above (on our nominal model) starting the value 
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of So randomly in the range [0, 2tt). The results obtained 
for 5d are shown in Fig. [9j The resolution extracted for 
So is 7.5°, which is less than the uncertainty on current 
measurements. Fitting with So free has no significant 
effect on our results for rs(£'A"*(892)), Sb{DK*{%Q2)) 
and 7 (for all model values of Sb)- 



80 



\i = -157.3° 




-180 



-160 



-140 



FIG. 9: Distribution of fitted results for Sd obtained in the 
amplitude analysis when fitting with 5d ss a, free parameter. 
The mean and width obtained by fitting the distribution to 
a Gaussian lineshape (solid line) are shown on the plot. The 
generated value of So is —158°. 



VI. EXPERIMENTAL EFFECTS 

The study reported in the previous section neglected a 
number of experimental effects which will surely impact 
the sensitivity to 7 in any real- world experiment. The 
effect most likely to have the largest impact on the res- 
olution is the presence of background events. To study 
how this might affect the sensitivity to 7, we generate 
background events by sampling from a flat Dalitz distri- 
bution. For real data, there will almost certainly be some 
features present in the shape of the background; how- 
ever, for our purposes, a flat background model should 
be sufficient. We take the expected background yields in 
each of the six flnal states to be equivalent. This should 
be a reasonable approximation provided the main source 
of background events is from combinatorics (a plausible 
assumption for LHCb) . The ratio of background to sig- 
nal, B/S, is defined here as the ratio of the number of 
expected background events to the number of expected 
signal events for the final state with the smallest expected 
signal yield. We generate background events according 
to this prescription for the ratios B/S = 1,2, 10, 50 and 
100, and fit the data using the same procedure as in the 
previous section. 

For B/S > 10, ambiguities in the solution begin to 
appear (see Fig. 10 1, just as they did using the quasi- 
two-body approach without background. About 2% of 



pseudo-experiments with B/S ^ 10 find their best likeli- 
hood in an ambiguous solution. This number steadily in- 
creases as the background increases, reaching about 25% 



for B/S = 100. This can be understood since as the 
background increases, it becomes more difficult to extract 
cleanly the amplitudes with smaller contributions. The 
amplitude with the strongest contribution in our model is 
DA'* (892). Thus, in the presence of a large background 
contamination, the amplitude analysis essentially reduces 
to the quasi-two-body approach. This results in the am- 
biguities in the solution found in the quasi-two-body ap- 
proach (without background) appearing in the amplitude 
analysis if the background yields are large. 

Our studies show that the background level should be 
kept to B/S < 20 in order to take full advantage of the 
benefits of the amplitude analysis. This appears achiev- 
able, even in a hadronic environment. For B/S = 20 
the ambiguities in the solution still appear at the few 
percent level; however, one can study the likelihood con- 
tours to determine whether the solution for any given 
experiment may have ambiguities. The results obtained 
from an experiment with such a background would need 
to take this into account. Note also that these results 
are for a data sample roughly equivalent to one nominal 
year's data taking at LHCb, and that effects due to am- 
biguities would be expected to be ameliorated with larger 
data samples. 

The results obtained for rpiDK* (892)), 
Sb{DK*{892)) 



XII 



The 



and 7 are given in Tab. 
ambiguities in the solution discussed above have been 
ignored in these results, i.e. the means and widths for 
B/S > 10 are obtained by fitting the region around 
the true solution to a Gaussian lineshape. The effects 
of the background on the resolution arc minimal. The 
resolution of 7 for B/S = 10 is only 1° worse than with 
no background. Even for B/S = 100 the resolution 
is only about 60% worse. Thus, the most important 
impact of the presence of background events is their 
ability to lessen the power the amplitude analysis has 
to break the ambiguities found in the quasi-two-body 
analysis. 



TABLE XII: Results for rB(DA'*(892)), SgiDK* (892)) and 
7 obtained from the Dalitz-plot analysis, with different levels 
of background. These results are for rs ~ 0.4, Sb ~ 0° and 
7 = 60°. 
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58.8 
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Another experimental effect which must be taken into 
account in a real experimental analysis is the variation 
of reconstruction efficiency across the Dalitz plot. Our 
study has assumed that the efficiency is flat, but in reality 
one would expect the probability to reconstruct success- 
fully a decay to be lower at the edges of phase space near 
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FIG. 10: Distributions of fitted values of 7 in tlie Dalitz plot approach, for B/S — 10, 50, 100. These results are obtained with 
7 = 60°, (5b = 0° and rs = 0.4. 



the corners of the Dahtz plot. If the shape of the effi- 
ciency function is known, this can be taken into account 
in the analysis. However, systematic effects arise since 
the efficiency is typically measured using Monte Carlo 
simulations, which inevitably will not give exactly the 
same behaviour as the data. 

It is impossible to give a quantitative estimate of how 
large an effect this may be. Nonetheless, there is a strong 
qualitative reason to believe that it will not be major 
problem. The key point is that the our amplitude anal- 
ysis extracts 7 from the difference between DKn Dalitz- 
plot distributions with the D meson reconstructed in dif- 
ferent decay modes {D Kir and D CP eigenstates) . 
Data/MC differences can be expected to cancel in the ra- 
tio of Dalitz-plot distributions, unless there is a momen- 
tum dependence in the ratio of efficiencies to reconstruct 
the D meson in the different final states. Such an effect 
should be possible to study using control samples. 

VII. SUMMARY 

We have presented a feasibility study of an extension 
to the recently proposed method to extract the Unitarity 



Triangle angle 7 from amplitude analysis of B ^ DKtt 
Dalitz plots. The analysis includes the cases where the 
neutral D meson is reconstructed in CP-even eigen- 
states as well as in CKM-favoured and CKM-suppressed 
hadronic decays. Compared to the previously proposed 
quasi-two-body analysis, the amplitude analysis provides 
(i) at least 50% better sensitivity to 7, (ii) resolution of 
ambiguous solutions, (iii) much reduced dependence of 
the sensitivity on the strong phase Sb, and (iv) the possi- 
bility to determine the poorly known parameter Sd- The 
analysis appears to be relatively robust against mismod- 
elling of the Dalitz plot, and performs well even when rel- 
atively large backgrounds are present. We conclude that 
this method appears to be a highly attractive addition to 
the family of methods that can be used to determine 7. 
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